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Addressed is the question of whether a natural mechanism exists to resolve the 
strong CP problem. The analogous issue for the two-dimensional CP N—1 models 
is analyzed using computer simulations. 



1 Introduction 

Non-perturbative effects can introduce CP violation into the strong interac 
tions. Since no such violations of this discrete symmetry have been observe 
an avoidance mechanism must be operative. Although many proposalaxl 
have been suggested to resolve this issue, the strong CP problem remains 
one of the outstanding low-energy mysteries of the standard model of particle 
physics. 

A four-dimensional Yang-Mills theory has instantons. They represent tun- 
nelling between different classical n-vacua, where n is an integer character- 
izing a topological winding number. The true quantum vacua ace believed 
to be linear combinations of the n-vacua weighted by exp(in9).au Alterna- 
tively, on can add to the standard Yang-Mills action the term 9Q, where 
Q = 5 2 /(32tt 2 ) / d 4 xF^F£"(x) is the total instanton number. The result- 
ing theory explicitly breaks parity, time reversal invariance and CP unless 
equals or n. When 9 equals ir, CP is believed to be spontaneously broken. 

When quarks are included in the analysis, the question of CP violation 
is more complicated. An overall phase in the quark mass matrix M. can be 
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eliminated through an axial /7a(1) transformation, but this produces an ex- 
tra contribution to 9 of Arg detM.. The relevant parameter is thus 9 e ff — 
9 + Arg detAi. In QCD, strong CP violation arises unless 9 e ff = 0. From 
measuremenla^on the neutron's electric dipole moment, a bound of e ff < 10~ 9 
is obtainedfrQ The strong CP problem can be stated as follows. Why is 
8 e ff < 10~ 9 when 9 e ff is a expected to be a parameter of order unity? In 
other words, 9 must be tuned to cancel the contribution to Arg detM to one 
part in a billion. If this is adjusted by hand at the tree level, loop effects ruin 
the cancellation. Various proposals have been made to resolve this problem, of 
which the most elegant is the Pecci-Quinn mechanismB It automatically sets 
6 e ff to zero. 

A number of years ago, I raised the question of whether the strong CP 
problem in pure Yang-Mills theories could be solved naturally.Q One should 
realize that it is strong-coupling dynamics that is most relevant to the issue 
of 9 dependence. Why is this so? Long-distance physics always determines 
vacuum structure. For example, local fluctuations in magnetic moments do 
not determine whether an iron bar magnetizes. A microscopic region may look 
disorderd but the overall system can be ordered, or vice-versa. Since 0-vacua 
are at the heart of the CP problem, the existence or non-existence of long- 
distance correlations is crucial. Due to our current limited understanding of 
long-distance physics and strong-coupling dynamics in Yang-Mills theory, per- 
haps some subtle effect concerning the strong CP problem has been overlooked. 

There exists a toy model in which the strong CP problem is dynamically 
resolved. It is the (2+l)-dimcnsional Georgi-Glashow model that was analyzed 
by Polyakov. It exhibits a natural relaxation mechanism by producing its own 
axion-like field to eliminate the 9 dependence. 

The model has instantons: They are the 't Hooft-Polyakov monopoles. 
Polyakov showed that these monopoles lead to confinement and produce a 
mass for the unbroken vector boson. One can repeat Polyakov's analysis^ in 
the presence of a 9 term. To the action, one adds 9Q, where Q is the net 
monopole number. The dilute gas approximation gives 

J dxd^^d^x) + 2\cos((f>{x)/^/JI+9) (1) 

for the partition function, where A and /i are parameters. Here 0, a field 
that is dynamically generated, reproduces the Coulombic interactions between 
monopoles. Since the shift <j> — > <j> — ^/JI9, eliminates the 9 dependence in Zg, 
the (2 + l)-dimensional analog of the strong CP problem is resolved. 

There is a physical reason why the 9 dependence goes away. Three- 
dimensional Coulombic gas systems are overall neutral. If Q = 0, then the 
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term #x(monopole number) has no effect when added to the langragian. 

But why is the three-dimensional Coulombic gas system neutral? The 
answer is due to the attractive 1/r potential between monopoles and anti- 
monopoles. The force between the monopoles and the anti-monopole is suffi- 
ciently strong so as to render the gas neutral. 

Once one realizes this, it is easy to establish criteria for when the 9 de- 
pendence in a theory naturally goes away.El Let V(x) be the anti-instanton- 
instanton interaction potential in a d-dimensional system computed in the ab- 
sence of any other instantons and let R be the distance between the instanton 
and the anti-instanton. 

If V(R) ~ l/R n with n < d then there is no 9 dependence . (2) 

For four-dimensional Yang-Mills theories, there is no strong CP problem if the 
instanton is attracted to the anti-instanton through a potential stronger that 
1/R 4 for R large. Because of our ignorance of the long-distance behavior of 
Yang-Mills theories, it is not known whether the criterion in Eq. (0) is satisfied. 



2 Computer Results of Schierholz et al for CP Models 

During the past few years, G. Schierholz and co-workers have performed some 
intriguing Monte Carlo, studies for the two-dimensional CP 1 ^^ 1 models in the 
presence of a 9 term.PH 12 ! These models have several features in common with 
d = 4 Yang-Mills theories: They have confinement and instantons. 

The free-energy difference per unit volume f(9) is determined from 

exp{-Vf(9)) = Z(6)/Z(p) , (3) 

where V is the volume of the system and Z(9) is the partition function with 
the two-dimensional 9 term of 9/(2ir) J d 2 xF$i. 

As a function of 9, the freepfinergy obtained from Monte Carlo simulations 
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exhibited a flattening behaviorliiliij that is reproduced well by 

U((3) 9 2 9 <9 C 
/_ \c(/3) 9>9 C ■ W 

In other words, / rises quadratically with 9, turns over and then displays no 
dependence on 9 beyond a critical value 9 C . For two-dimensional confining 
systems, the string-tension a for a particle of charge e can be computed from 

a(e,9) = f(9 + 2ne)-f(9) . (5) 
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Thus the behavior in Eq. (Q) implies that confinement is lost for 9 > 9 C for 
sufficiently small charges e. Furthermore, the Monte Carlo simulations sug- 
gested that 9 c (j3) goes to zero as the continuum limit is taken. Taken together, 
the above statements imply that a confining continuum CP 1 ^^ 1 theory would 
have to have 9 adjusted to zero. The conclusion is that the requirement of 
confinement necessitates a solution .of the two-dimensional analog strong CP 
problem in the CP N ~ 1 models.BB 

By the way, a flattening behavior in f{9) was also seen for the d = 4 SU(2l 
lattice gauge theory in subsequent simulations performed by G. Schicrholz.E3 
In d = 4, however, this does not imply the lost of confinement so that the 
intepretation of this result is not so clear. 

Inspired by Schierholz's work, Jan Plcfka and I decided to perform our 
own studies. We were particularly intereseted in whether the "deconfinement 
effect" could be related to a natural relaxation method. 



3 General Discussion of Lattice Measurements of the Free-Energy 
for a System with a 9 Term 

It is not easy to simulate a system with a 9 term. This is because one usually 
uses the 9 — vacuum to analyzed the 9 ^ vacua. If "barrier penetration" 
effects are strong, poor results can ensue. 

The free-energy difference f(6) can be computed from 

e3tp(-V/(0))=$^P(Q)exp(i0Q) . (6) 



Here, P(Q) is the probability that the system has topological charge Q. Let 
Nmc{Q) be the number of times that Monte Carlo configurations with topo- 
logical charge Q are generated. Then 

p / n , _ N M c(Q) (7 s 

Pmc{Q) = E^WQO (7) 

provides a Monte Carlo estimate Pmc of P(Q)- Using Eq. ([]]), a Monte Carlo 
value /jifc(S) of f(9) can be obtained from 

exp(-Vf M c(9)) = ^PMc(Q)eMWQ) ■ (8) 
Q 

This method of computing f(9) can lead to an anomalous flattening effect. 
Suppose there is a dominant statistical error in the Q = sector and that the 
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probability of having^ Q — configuration is over estimated: Pmc(0) > P(0). 
Then one can showtLa that 

fMC ~\f(6 b ) e>e b > ^ 

where 5P(0) = P M c{0) - P(0) and 

f(9 b )^±\\og\SP(0)\\ . (10) 
Here are the key points (see reference [14]): 

(1) When the volume is sufficiently big, a limiting 9 b arises. For 9 > 9 b , 
f(6) cannot be reliably measured. 

(2) If a flattening behavior in f(9) is seen, one should be suspicious of the 
result. One needs to check whether |<5P(0)| is bigger than it should be. 

(3) If a large- simulation shows a flattening effect for f(9), but a smaller- 
V simulation does not, then one should probably trust the smaller-]/ result. 

In refs.[ll,12], no flattening effect was seen in simulations with small vol- 
umes. Point (3) says that one should trust the smaller volume studies. If true, 
the CP N ^ X models do not undergo a deconfining phase transition. 

4 Monte Carlo Results for the Exactly Solvable d — 2 U(l) Gauge 
Theory 

One way to decide the question of anomalous flattening behavior is to study 
an exactly solvable system. The d — 2 U(l) gauge theory is such a case. The 
action is given by 

s u w=pJ2(u p + u;) , (ii) 
p 

where U p is the product of the U(l) link phases around a plaquette p and (3 is 
the inverse coupling. 

The lattice analog of the continuum (9-term action i9/2n J cPxFqi is 

term log (Up) . (12) 

P 

The exact analytic result for the partition function governed by the sum of the 



action in Eqs. (11) and (121) has been computed in ref.[15]. 

When Jan Plefka and I performed Monte Carlo simulations of this system, 
we sometimes found a flattening behavior in the free energy that disagreed 
with exact analytic results. lii Figure 1 illustrates this. 
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Figure 1: The Free Energy of the U(l) Gauge Model Versus 8. The solid line is the exact 
result. 

Noteworthy features are: 

(1) A barrier 8b arises. 

(2) A short calculation verifies that 9b is the value expected from Eq. ([To|). 

(3) The flattening effect is due to P MC (0) > f(0). 

ft turns out that Ob moves to smaller values as V increases. For a 7 x 7 
lattice, the free energy agreed with exact analytic results for all f(6). 

5 Monte Carlo Results for a CP*' 1 Model 

We also decided to simulate the CP 3 model for a lattice action involving an 
auxiliary U(l) gauge field. The advantage -af-using this latticization is that 
analytic strong coupling series are available JiaO The action is given by 

S = P N ^2( z x- Zx+aU (x,x + A) + z x ■ z* +A U* (x,x + A)) , (13) 

x,A 
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Figure 2: The CP 3 Free Energy Versus 9 at /3 = 0.2. The solid line is the strong coupling 
series. 



where the complex scalar fields z % x satisfy Yli=i z xi z x = 1 an d N — 4 for the 
CP 3 model. Here, the sum is over unit shifts A in the d positive directions, 
where d — 2. The field U (x, x + A) is a phase associated with the link between 
x and x + A: it is the U(l) auxiliary gauge field. Eq. ( [l3| ) and Sg of Eq. (12) 
are added to obtain the full action. This lattice action differs from the one 
used in refs. [11,12] but is equally good. 

When we performed Monte Carlo simulations, again, spurious flattening 
behavior in the free energy sometimes occurred. Figure 2 is an example. 



6 Conclusion 

Here are the main points: 

(1) Subtleties arise when simulating systems with a 9 term. The basic 
difficulty is that one is trying to simulate one vacuum (a 9 ^ vacuum) 
using another vacuum {9 = vacuum). For good Monte Carlo results, barrier 
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penetration needs to be suppressed. When the volume of the system is large, 
the barrier penetration becomes sufficiently severe that reliable results cannot 
be obtained, particularly for 9 near -n. 

(2) Because current simulation methods tend to have the biggest error in 
the measurement for the Q = probability, Monte Carlo results can generate 
an anomalous flattening behavior in the free-energy for 9 large. 

(3) Reliable results for the free-energy appear to be obtainable throughout 
the entire 9 region as long as the volume is keep sufficiently small. 

(4) The ideas and results of this presentation suggest that the CP 1 ^^ 1 
model with a 9 term does not undergo a deconfining phase transition. 

(5) In short, it is unlikely that the solution of the strong CP problem as 
suggested by previous numerical studies in the CP 3 model actually works. 
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